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Abstract 



In the following work we consider the possibility of explaining the observed 
baryon number asymmetry in the universe from simple baryon number violating 
modifications, involving massive scalar bosons, to the Standard Model. In these 
cases baryon number violation is mediated through a combination of Yukawa and 
scalar self-coupling interactions. Starting with a previously compiled catalogue 
of baryon-number violating extensions of the Standard Model, we identify the 
minimal subsets which can induce a B — L asymmetry and thus be immune to 
sphaleron washout. For each of these models, we identify the region of parameter 
space that leads to the production of a baryon number asymmetry of the correct 
order of magnitude. 
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1. Introduction and Motivation 



The excess of matter over antimatter observed in the present day universe 
can be explained by proposing that there existed suitable baryon number vi- 
olating physical processes which operated in the early stages of the universe. 
These physical processes led to the situation whereby the production of baryons 
was slightly greater than the production of antibaryons thereby producing the 
required baryon-antibaryon asymmetry. In fact the quark-antiquark asymmetry 
required in the early universe (t < 1CP 6 sec) to give rise to the present excess 
of matter over antmatter is very small: 

^^i~3xl0- 8 . (1) 
n q 

Thus the processes leading to the quark-antiquark asymmetry need only be 
minute to account for baryogenesis. 

The three basic requirements for baryogenesis, which were first compiled in 
[JD, are: (a) Baryon number violation: if baryon number was not violated then 
the present day baryon asymmetry could only be explained by an initial baryon 
number asymmetry; (b) C and CP violation: even with the existence of baryon 
number violation, a baryon asymmetry would not have arisen unless there was 
a preference for the production of matter or antimatter in the early universe; 
and (c) Non equilibrium conditions: if the universe were in a continual state of 
thermal equilibrium, the phase space densities of baryons and antibaryons would 
neccessarily be equal, owing to the fact that by CPT invariance the baryon and 
antibaryon masses are always equal. 

RefQ provides a catalogue of all of the simplest extensions of the Standard 
Model that explicitly violate baryon number. (With lepton number broken ex- 
plicitly through neutrino Majorana masses, all of the global U(l) symmetries of 
the Standard Model are broken and an understanding of charge quantisation re- 
sults H . The charge quantisation issue provides a strong theoretical motivation 
for constructing all of the simplest Standard Model extensions which explicitly 
break B. See RcfQ for further discussion.) In each of these models baryon 
number violation is mediated through scalar bosons, which are required to be 
massive so as to ensure that nucleon decay remain unobserved. In the following 
work we will be considering the possibility that the out of equilibrium decays of 
these massive scalar bosons in the early universe could account for the present 
day baryon number asymmetry. This work differs from similar calculations 
arising from GUT models Q in that most of the models we will consider break 
baryon number through a combination of Yukawa interactions and scalar self- 
couplings rather than through the baryon number violating Yukawa interactions 
associated with a single massive particle. 

Starting with the baryon number violating interactions proposed in we 
are left with the task of both demonstrating the existence of and calculating the 
magnitude of the CP violation arising in each of the baryon number violating 
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models. In each case there will be two opposing decay paths for our massive 
scalar each giving rise to products with different baryon numbers. We can 
analyse this type of system using the following procedure: 

If we assume for simplicity that our scalar boson X has the two quark lepton 
final states ql(B = —1/3) and qq(B — 2/3), then C and CP are violated if the 
branching ratio of the X to the ql final state (= r) is unequal to the branching 
ratio of the X to the ql final state (= f), i.e. r ^ f. We also know that CPT 
invariance requires that the total decay rates of X and X be equal. Therefore we 
can write the branching ratios of the corresponding decays X — > qq and X — > qq, 
as 1 — r and 1 — f respectively. From a symmetric initial state consisting of X 
and X, provided there are no further baryon number violating reactions, a net 
baryon asymmetry will exist after all of the X and X decay, with the mean net 
baryon number produced being given by, 

AS = ~r + ~(1 - r) + if - ~(1 - f) = -(r - r). (2) 

We can therefore evaluate the mean net baryon number produced by calculating 
the difference in the branching ratio for boson and antiboson decay for just one 
of the decay channels. 

We are hence essentially left with the task of calculating r — r, which is a 
measure of the CP violation in the system. CP violating effects come from con- 
sidering the interaction of the lowest order decay diagrams together with their 
loop corrections. It can be shown that the intermediate states in the loops must 
not only have CP violating complex couplings, but must also propagate on shell, 
thereby leading to complex Feynman amplitudes. Another important require- 
ment is that the loop corrections involve baryon number violating interactions 
(see H). In general we expect AB to be of order a N , where a characterises 
the coupling constants of the loop particles, and N is the number of loops in 
the lowest order diagram which interferes with the tree graph to give a non zero 
value for AB. 

2. The Model: 

The following modifications to the standard model were devised with the 
primary aim of obtaining complete electric charge quantisation from the gauge 
invariance of the Lagrangian. 

Complete charge quantisation through gauge invariance can be obtained 
provided there is only one unembedded U(l) invariance, and we obtain the 
correct charge quantisation provided the generator of this U(l) invariance is 
the standard weak hypercharge Y. If there is more than one unembedded U(l) 
invariance then the actual weak hypercharge of the theory can be chosen to be 
some linear combination of the standard model hypercharge and these additional 
symmetries of the theory. 
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As it stands the three-generation minimal standard model has five U(l) 
invariances, the standard weak hypercharge Y, baryon number B, and the three 
family lepton numbers L e , L^, and L T . These five U(l) invariances correspond 
to there being four classically undetermined electric charges. To remove these 
invariances we must construct extensions of the minimal standard model which 
break B and but leave Y exact. 

The simplest and most interesting way to break £7(1)^ is to introduce non- 
zero neutrino masses. This is most easily done by introducing right handed 
neutrinos into the model, where we choose that our right and left handed neu- 
trinos be related through Dirac and Majorana mass terms. 

This leaves us with just one undetermined electric charge, which can be 
taken to be the electric charge of the down quark Q(d), or equivalently the 
hypercharge of the down quark y d . Our four parameter uncertainty has thus 
been reduced to a two parameter uncertainty by this simple extension of the 
lepton sector. For more information see [Q. 

We are therefore left with the task of breaking this unwanted baryon symme- 
try without affecting the U(l)y hypercharge symmetry. This dual requirement 
is acheived by extending the standard model further to include either one or two 
exotic scalar multiplets together with a set of baryon number violating Yukawa 
and scalar-scalar self interactions. For further information on these models see 

§■ 

The Yukawa interactions and the SU(3)c®SU(2)l<3U(1)y representations, 
where the hypercharge Y is parameterised in terms of the unknown down quark 
hypercharge y d , is listed below: 
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Note that we have used the following notation for the standard model fermions 
and right handed neutrinos: 

f L ~ (1,2,-1), e R ~(l,l,-2), ^-(1,1,0), 

Ql ~ (3,2,1 + 2/d), u R ~ (3,1,2 + y d ), d B ~ (3,1,0,,). (4) 
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01.1 


= CT 3.1 


01.2 


_C 

- a 3.2 


04 


= a 7.1 



The fermion interactions, Ql(Ql) c , Ql{Ql) c , ur(ur) c , and dn(dR) c associated 
with the (73.2, (73.3, (76.1, and (77.1 scalars are flavour antisymmetric. 

The simplest baryon number breaking scalar interactions were obtained by 
noting that certain pairs of scalars in the above list have group properties which 
can be related by conjugation subject to the required constraint yd = —2/3. 
For these cases the two distinct scalars listed in Eq(^) can be considered as 
representing one particle with two sets of baryon number breaking Yukawa in- 
teractions. There are three such conjugate pairs which we list below: 

(3,1,2/3) 
(3,3,2/3) 

(3,1,-4/3). (5) 

Baryon number violation is also obtainable by considering the scalar-scalar 
interactions of the scalars listed above. By considering every possible scalar 
combination in Eq(^), two lists of possible charge quantising scalar potentials 
were compiled, corresponding to AB — 1 baryon number violating processes 
and AB = 2 baryon number violating processes respectively. The AB = 1 list 
is shown below: 

01.02 — ► 01.201.2020 

01.03 — > 01.103.1 + 01.101.101.103.1 + 01.103.103 .103.1 + 01. 103. 10 
— > 01.203.2 + 01.201.201.203.2 + 03.203.203.201.2 + 01.203.20 V 

(71,(75 — > (71.1(75(75 

— ► 01.201.2050° 

01.06 — ► 01.206.100 

01.07 — ► 01.207.10 C C 

02.03 — * 0203. 203. 20 C 

02.07 — * 0207.10 

03.04 — > 03.20400 

03.05 — > 03.1050 

— > 03.203.2050 

03.08 — > 03.2080 C C 

(74,(77 — > 0407.1 +04040407.1 + 07.107.107.104 + 0407.1 

05, P ~> 050505P 

05,05 0505050 

05.07 — ► 0507. 10 C 

06.08 — * 06.108 + 06.106.106.108 +08080806.1 +06.1080^0 (6) 

where represents the SM Higgs scalar <j> ~ (1,2,1), and p represents a new 
Higgs like scalar p ~ (8, 2, 1). Similarly the AB — 2 list of scalar potentials take 
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the form: 

103.1(71. 1CT3.1 

1C3.2C1.1cr3. 2 
1C3.3C1.11T3. 3 
1C3.4C1.1CT3. 4 
2C3.1C1.2C3.1 
2C3.2C1.2C3. 2 
2C3.3C1.2C3. 3 
2C3.4C1.2C3. 4 
1C3.1C7.2 
2C3.2C7.2 
3C3.3C7.2 
4C3.4C7.2 
C4CT7 — > C4CT7. 1C4CT7.1 

— » C4(T7.2C4(T7.2 
C6,C7 — » C6.2C7.1C7.1 

— ► C6.2C7.2C7. 2 . (7) 

The above scalar potential terms can be placed into groups consisting of 
quadratic, cubic and quartic terms. As shown in [0 each of these different 
groups of scalar interactions are of different phenomenological interest, with the 
higher order interactions (i.e. the quartics and cubics) having the least stringent 
experimental constraints (obtained from nucleon decay data) on m a , the mass 
of the scalar. Each of these subgroups will also result in different outcomes as 
far as explaining baryogenesis is concerned. In the following work we will be 
determining which of the above interactions are able to account for baryogene- 
sis. This will basically involve the calculation of r — f, the CP violation arising 
from each model and the subsequent use of standard cosmological arguments. 



ci, c 3 



C3, c 7 



ci. 
ci. 
ci. 
ci. 

Cl. 

01. 
ci. 
ci. 

C3. 
C3. 
C3. 
C3. 



3. Calculation of Baryon Production: 

As previously mentioned baryogenesis is only possible provided we have CP 
violation incorporated in such a way as to induce different partial decay rates 
for baryon number violating decays of particles and antiparticles. To obtain 
this CP violation we require the introduction of at least two baryon number 
violating scalar self interactions. For example for the ci. 105(75 class of scalar 
potentials which we will consider later, we are required to introduce two scalar 
interactions, ci. 1(75(75 and ai.\a' 5 a' 5 , where the scalars (75 and a' 5 have identi- 
cal group properties but different masses and couplings. There are two reasons 



G 



behind this requirement; the first being that if the exchanged scalar (s) Y in 
the loop corrections have identical masses and properties to the decaying par- 
ticle^) X, then the contribution to baryogenesis arising from the decay of X 
with the exchange of Y will be cancelled by an equal and opposite contribution 
made by the decay of Y with the exchange of X; the second reason arises from 
the requirement that there be an imaginary component to the product of cou- 
pling constants associated with the CP violating tree and loop graphs, which is 
impossible if each Yukawa coupling appears with its conjugate. 

The above list of baryon number violating models can be narrowed down 
if we consider the effects of sphaleron or other forms of damping processes on 
any baryon number asymmetry produced by each model || . The rapid B + L 
violating sphaleron transitions that are still occuring after the decays of our 
o-particles will quickly erase any B + L asymmetry. If we take this factor into 
account we can rule out a number of the models listed in Eqs(||||Q), by keeping 
only those that produce a B — L asymmetry (these being immune from sphaleron 
washout) we are left with the following much shortened list of possible baryon 
asymmetry generating models: 

Ol,0 5 — » O1.1O5O5 
0-2,(77 — > O2O7.10 
03,05 — > O3.1O50 
05,07 — > O 5 O7.i0 C 

(8) 

Thus if we allow for damping, we are left with just four interactions, one cubic, 
and three Higgs containing cubic interactions. 

It should also be noted that the version of the SM which we have used 
includes Majorana masses for neutrinos, which can in themselves allow for the 
production of a baryon asymmetry or the damping of a preexisting baryon 
asymmetry. If M > m a , where M is a Majorana neutrino mass, then out- 
of-equilibrium decays of vr may produce a AL prior to the decays of the o 
bosons. However, this asymmetry will be erased, and in particular it will not 
be reprocessed into a AB through B + L violating sphaleron processes. This 
is because the B — L violating o interactions (which we assume still occur 
rapidly) will combine with sphaleron processes to force both AL and AB to 
vanish. If M < m a , then any AB produced by o decays will be erased by the 
combination of rapid L violating and B + L violating processes. Subsequent 
out-of-equilibrium decays of vr may produce a AL which gets reprocessed into 
a AB by sphaleron effects, given that the o-induced B — L violating processes 
have already switched off. Although this is an interesting scenario (and has been 
studied in other contexts 0), it is not the one we wish to consider in this paper. 
We will therefore require that M > m a , so that only o-decays contribute to the 
AB that survives to the present day. (We also need to assume that no baryon 
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asymmetry is produced during the electroweak phase transition. This, however, 
is assured provided that we assume the existence of only a single electroweak 
Higgs doublet <f>, see for example RefJ^] for a review.) 

Although any baryon number asymmetry arising from the conjugate pair 
and quadratic interactions will be damped away, we consider these systems first 
as an instructive exercise in how to do the calcuations. The realistic cases will 
then follow by simple extension. 

The conjugate pair interactions are the simplest class of interaction included 
in our catalogue. The calculation of the baryogenesis arising from these conju- 
gate pair interactions is directly analogous to that obtained by ||| using a GUT 
model. CP violation arises from considering the tree and one loop corrections 
for the scalar decay 0. For the ai.i — of i system, our Yukawa terms have the 
following form, 



£ = \i{e R ) c <ru R + \ 3 u R a(d R ) c 

+ X' 1 (e R ) c a'u R + X 3 u R a'(d R ) c + h.c. 



(9) 



where a ~ a' ~ (3, 1, 2/3)_under SU(S) C ® SU{2) L <g> U(1) Y . The reasons for 
introducing two copies of (3, 1, 2/3) are discussed in the first paragraph of Sec 2 
above. This leads to the following expression for the baryon number production: 



4ImTr(X\X[X 3 X' 3 ') 
Tr(X\\i) + Tr(X\X 3 ) 



Iml(m a i /m<r)- 



(10) 



The function I(m a : /m a ) represents the Feynman integral for the exchange of 
the scalar a' in the decay of the scalar cr, and it has the form 



Iml(m a i /m a ) 



1 

16tt 



1 - 



In 1 



(11) 



If we assume that the mass of a is much larger than the mass of a' then we 
can ignore the contribution made by the decay of the a' scalar and Eq(|l^) is 
essentially the complete contribution to baryogenesis. 

In the quadratic models baryon number is violated through a scalar-scalar 
interaction. The analysis of the bilinear interaction coupling two different scalars 
is similar to that associated with the conjugate pair case, and we must again 
introduce two sets of interactions to obtain non zero CP violation. In our 
sample calculation we will use the <J\.\o 3 .\ and a[ ^3 1 bilinears together with 
the following Yukawa interactions: 



£ = Xi(e R ) c aiU R + \ 3 (d R ) c a3,u R 
+ X' 1 (e R ) c a' 1 u R + X 3 (d R ) c a' 3 u R 



(12) 



h.c. 



The CP violating tree and one loop corrections for this model are shown in 
Fig(l). In this case we have a four vertex loop correction rather than the 
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three vertex loop correction we had for the conjugate pairs. The CP violating 
Feynman amplitude corresponding to o\ decay takes the form: 

M{o\ — > u + e) = Ai + ■^-^—^\[\t\' 3 I(m l7 > /m (71 ,m CT '/m CT1 ), (13) 
ml — mi 1 3 

where /z 2 is the coupling constant associated with the bilinear (J, a±, 1(73,1, and /i' 2 
is the coupling constant associated with the bilinear p' 2 <j[ i<y' 31 , where p and p! 
have dimensions of mass. The factor J(rn CT ' /m (T1 , m CT / /to ct1 ) again represents the 
contribution of the Feynman integral around the loop; in this case two scalars 
connected by a mass insertion are exchanged. By taking the difference between 
absolute value of Eq(fi"3"|) squared and the corresponding squared amplitude of 
the antiparticle process, we arrive at the following expression for the baryon 
number production: 

4 ImTrUp?)^^^) 

(14) 

The value of I(m a > /m ai , to ct ^ /m ai ) in this case is found to have the form: 



Iml(p 1 ,p 3 ) = on _ 2 _ 2 (-7pj + lp\ + p\ ln[l + p\ 



1 1 

32n p\ — p\ 

7pf ln[l + P r 2 ] - Pi ln[l + Ps 2 ] - 7pt ln[l + ft 2 ]) (15) 



where, 



= (16) 



Note the similarity between the two expressions for the imaginary components 
of the phases for the conjugate pair and quadratic cases, Eq@ and Eq©, 
particularly in the limit m a > /m a — » oo. We assume that m CT3 > m CTl > m CT j 
which as can be seen from Eqs(|l5|) and (jL4|) will allow us to ignore the contri- 
bution the decays of a' 13 and 03 will make to baryogenesis. Thus Eq(|l4|) can 
be regarded as constituting the entire contribution to baryogenesis. 

It should be noted that we can also obtain baryogenesis in the special case 
where we have just three new exotic scalars. For example for <j\ decay we can 
set (T 3 = (T3 and still obtain the required imaginary components to the Yukawa 
couplings and the Feynman integral. 

We now turn to realistic models based on Eq(||). In addition to not having 
a sphaleron washout problem, these models are also of greater phenomenolog- 
ical interest because the nucleon-decay bounds on m a are much weaker than 
those for the unrealistic conjugate pair and quadratic toy models just consid- 
ered {m a > 10 6 GeV for the a%. 1(75(75 cubic compared with vn a > 10 15 GeV for 
the conjugate pair and quartic interactions 0). 
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From a diagrammatic and calculational point of view these models differ 
from the above because baryon number violation arises from a cubic or trilincar 
term in the Higgs potential. 

Consider first the <T\.\ — 05 system. We will initially introduce two copies of 
both the CTi.i and 05 scalars; we will require at least three of these four scalars 
to obtain the required CP violation. We will consider a simplified set of Yukawa 
interactions associated with the participating scalars as shown below: 

C = \i{e R ) c <Jiu R + X 5 v L a 5a d R + X 5 e L cr 5b d R (17) 
+ A'i (e fl ,) c 0i ur + X' 5 D L a' 5a d R + X' 5 e L cr' 5b d R + h.c. 

where 05a and 5 {, designate the two SU(2) components of 05. The above set 
of Yukawa interactions are simplified in that we have considered just a few 
of the Yukawa interactions, see EqQ, which may from group considerations 
be associated with the exotic scalars. The CP violating tree and loop corrected 
diagrams for the decay of the 05 multiplet are shown in Fig(2). For these decays 
the lowest order loop corrections involve both the 010505 and 010505 cubic 
interactions. (In the interests of simplicity we have ignored the contribution 
the 0i050 5 interaction will make to baryogenesis.) Thus we require only one 
version of the 01 scalar to obtain baryogenesis from the decay of the 05 scalar, 
conversely we would require two versions of the 01 scalars and one 05 scalar 
to obtain baryogenesis from the decay of the 01 scalar. The loop corrected 
diagrams shown in Fig(2) are complicated and as such won't be evaluated. 
This omission is justified because we are interested in an order of magnitude 
calculation only. In any case, it is the coupling constants which will play the 
major part in determining the numerical value of the baryogenesis arising from 
these cubic interactions. Based on the results obtained for the conjugate pairs 
and the quadratic interactions it is reasonable to assume that the imaginary 
part of the Feynman integral will be of a similar form and hence give similar 
numerical values to those calculated in Eqs(ll) and (|l^); we must however allow 



for an extra l/(2ir) 2 suppression factor for each additional loop order. For the 
diagrams given in Fig(2) the imaginary component of our loop integral will 
thus be suppressed by an extra factor of l/(27r) 2 in comparison to the one loop 
integrals considered in Eqs(|ll|) and (|l5|). The CP violating amplitudes for the 
decays shown in Fig(2) are given by: 

M(a 5 ->d + /) = A 5 + fifi" 1 X' 5 \lX' 5 I (m CTl ,m ff5 ,m^). (18) 

where /i represents the coupling constant associated with the cubic interac- 
tion (again with units of mass), and I(m ai , m CT5 , m a i ) represents the Feyn- 
man integral taken around the loops, which in this case will not be evalu- 
ated. From Eq(|l8|) we obtain the following expression for the total baryon 
number production arising from the decays of the 05a and 05b scalars, where 
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m <?5a — m <?5b I 

ktv \t 



Note that we have included an extra factor of two in the above expression to 
allow for the fact that we have two equal (as m aha ~ m a 5b ) contributions to 
the baryon number production. If we assume that m CT5 > ra„; then we can 
ignore the contribution the decay of a' 5 will make to baryogenesis. Thus Eq([l9|) 
is in effect the complete contribution the system given in Eq(|l7|) will make to 
baryogenesis. 

The phenomenologically less interesting Higgs containing trilinears (m a > 
10 11 GeV j2j), will lead to baryogenesis via CP violating loop diagrams such as 
that shown in Fig(3) for the a 3 ,ia^(j) cubic, where we have used the following 
Yukawa interactions, 

£ = h(uR) c cr 3 d R + X 5 v L a 5a d R (20) 
+ A 3 (u fl ) c cr 3 d fl + X' 5 P L a' 5a d R + h.c. 

Note the similarity between the diagrams in Fig(2) and Fig(3). For er,; decay 
resulting from the interactions \x<j 3 a^4> and /i'a 3 a' 5 (f> we obtain an expression for 
r — f of the form, 

r -r = —~ — — Iml(m a5 ,m a3 ,m a5l , m^^rn^). (21) 

mi, Tr(XlX 5 ) 

If 0-5 is significantly more massive than 03 , a' 3 and a' 5 then Eq(|2l]) is essentially 
the complete contribution to baryogenesis. In this case our expression for the 
Feynman loop integral will involve the Higgs scalar with m a . 

Note that in the special case of 173 = a' 3 we can still get baryon number 
violation from 05 decays, and vice-versa for the decay of the 173 scalar. Thus 
as in the trilinear case the simplest possible baryon number violating system 
requires the introduction of two Higgs containing cubic interactions and three 
exotic scalars. 



4. Numerical estimates 

From our expressions for the baryon number creating r — f obtained in the 
previous section, we will attempt to numerically determine the likelihood of each 
model accounting for baryogenesis. 

We know from |J that the decay of super heavy bosons and antibosons at 
temperatures sufficiently well below their masses will produce a cosmic baryon 
entropy ratio of the form; 

— - = 0.28— AB, 22 

s N 
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where k is Boltzmann's constant, Nx /N is the ratio of the number of boson he- 
licity states to the number of light particle helicity states (Nx/N ks 10~ 2 ) and 
AB represents the baryon number generation arising from the decay. Astro- 
nomical observations put fcns/s « 10 -10 to 10 -8 , therefore AB must lie within 
the range 10~ 8 to 10 -6 , to account for the presently observed baryogenesis. 

Our expression for the baryon number created by the cubic interactions 
Eq(|l9|) can be expressed as 

AB = ^ ^ e\ 2 lml(m a , m a , m a > ). (23) 

Similarly our expression for the baryogenesis arising from the Higgs containing 
cubics, see Eq(^l]), can be expressed as 

AB = — — eA lml(m a . , m a3 , m a > , m a > , rrid,) (24) 

where A represents a typical value for the magnitude of the unknown Yukawa 
couplings in Eqs(|l^) and (^l]), and e is a phase angle characterising the aver- 
age strength of the CP violation associated with our scalar-scalar, and Yukawa 
interactions. 

In our analysis we will assume that the scalar self interacting and Yukawa 
couplings are approximately equal in magnitude. Using this assumption we 
will set out to find the region of parameter space within which these coupling 
constants must fall in order to account for baryogenesis. 

The value of the imaginary component of the Feynman integral for the above 
two cubic interactions is estimated from the corresponding expressions obtained 
from the conjugate pairs and quadratic interactions. After allowing for the 
extra l/(27r) 2 integration factor, it is thus assumed to fall within the range 
(27r)~ 2 (10~ 3 — 10~ 2 ). The value for the lower bound arises from the equal mass 
m a = °n\ a i and pi — 1 values of Eqs( |ll|) , and (^|) respectively, and the upper 
bound arises from the consideration of the m a 3> ra a i limit of Eq([H]) . 

By taking all of these assumptions on board (for both classes of cubic in- 
teraction), it is found that the coupling constants A and [i must lie within the 
following range to account for baryogenesis, 

Ae l/4 ~ M e l/4 « 1Q -1 _ 1Q -2 ( 25) 

If we assume for example that 1 > e > 10 -3 rads, then the couplings A and 
\l±\/m a must be of order 1 — 10~ 2 to account for baryogenesis. 



4. Conclusion: 
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We have thus demonstated that baryogenesis can be induced by each of the 
model types listed in Eq(||). 

In the numerical analysis of the various interactions listed in Eq(Q), it 
was found that both classes of cubic interactions can successfully account for 
baryogenesis provided the coupling constants fall within the range, Ae 1 / 4 = 
/le 1 / 4 /?/^ « 10 _1 — 10~ 2 . Given that the known Higgs Yukawa couplings range 
up to order unity for the top quark, it is not at all unreasonable that the un- 
known Yukawa's may fall within this parameter space. 

We have thus demonstrated that at least some of the baryon number vi- 
olating models introduced in 0], with the explicit aim of obtaining complete 
charge quantisation from gauge invariance will also allow for baryogenesis. It 
should however be noted that whilst we only need one of the baryon number 
violating scalar interactions listed in Eqs(||,||j^) (and consequently two exotic 
scalar multiplcts), to obtain charge quantisation, we require the introduction of 
at least two scalar-scalar interactions and at least three exotic scalars to account 
for baryogenesis. Thus the minimal requirement for charge quantisation is not 
sufficient to also give us baryogenesis. 
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Figure Captions: 



Figure 1: 

CP violating tree and one loop corrected diagrams resulting from the decay 
of u\,\ in the 0-1.10-3.1 quadratic model. 

Figure 2: 

CP violating tree and loop corrected diagrams resulting from the decay of 0-5 in 
the 0-1.1CT5CT5 cubic model. 

Figure 3: 

CP violating tree and loop corrected diagrams resulting from the decay of 0-5 in 
the o- 3 .iCT 5 cubic model. 
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figure 2. 
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figure 3. 
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